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Abstract. The absorption of ultraviolet light creates excitations in DNA, 
which subsequently start moving in the helix. Their fate is important for an 
understanding of photo damage, and is determined by the interplay of electronic 
couplings between bases and the structure of the DNA environment. We model 
the effect of dynamical fluctuations in the environment and study correlation, 
which is present when multiple base pairs interact with the same mode in the 
environment. We find that the correlations strongly affect the exciton dynamics, 
and show how they are observed in the decay of the anisotropy as a function of a 
coherence and a population time in a non-linear optical experiment. 



1. Introduction 

The bases that form the DNA molecule absorb ultraviolet light, leading to a surplus 
of energy that can cause damage. Fortunately for all living organisms that depend 
on this molecule, the excess energy is dissipated rapidly. The existence of delocalized 
states has been suggested to play an important role in the dissipation mechanism 
[I]. Evidence for delocalized excitations was found using femtosecond time-resolved 
spectroscopy [2j [3], as well as from combined density functional and molecular 
dynamics simulations [H [5] . 

Spectroscopic studies reveal a shift of 5000-7000 cm^^in the energy of an exciton 
in DNA after excitation. This shift is observed as a Stokes shift, the difference between 
the maxima of the absorption and steady state fluorescence spectra [5], as well as by 
a time-resolved shift in the excited state absorption [3]. In the latter experiment, the 
wavelength of the excited state absorption is found to shift from 380 to 330 nm, which, 
in the case of localized final states, is equal to the shift in the energy of the created 
excitation. 

This observed energy shift could be explained in a band model as follows. Because 
the transition dipoles in stacked bases make only a small angle (36 degrees in B- 
DNA) with each other, and positive values have been reported for the couplings, the 
DNA stack forms an H-aggregate. In this case, absorption is mainly at the top of 
the excitonic band, while emission occurs from lower-lying states. This qualitatively 
explains the observed shift. Quantitatively, however, the maximum Stokes shift in 
such a model is given by the exciton band width. As a rough estimate, which does not 
include the long-range interactions, this is equal to four times the coupling strengths. 
Most reported values of the couplings are smaller than 300 cm~^. Even when the 
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long-range interactions are included, significantly larger couplings would be needed to 
explain the observed Stokes shift. More naturally, however, the shift can be explained 
from an interaction with the environment, which increases the average bandwidth. 

Such an interaction has been included by considering the effect of a static 
disordered environment on the dynamics of excitons and charge transfer states [3 H] . 
In the case of charge transport in DNA, which is often interpreted in band models [9], 
the effects of disorder have been studied extensively. The effect of a static environment 
on the excitation energies [10] , as well as static variations in the structure of the DNA 
[TT] have been included in calculations of the conductance. Studies have highlighted 
the role of correlated fluctuations, present when multiple bases in DNA interact with 
the same mode in the environment, and their effect on the electronic properties p^ll3) . 

In the static description, the environment does not change on the time scale of 
the experiment. Its state is, however, different for each helix in the ensemble probed 
in a measurement. In turn, the system properties such as the transition energies of 
each base, or the couplings between bases, vary from helix to helix. Their values for 
an individual system are often called realizations of the disorder, and the properties 
of the ensemble can be found by averaging over all possible realizations, weighted by 
their respective probabilities of occurrence. In the static situation, the eigenstates of 
the Hamiltonian are well-defined, and can be used as the first step in understanding 
the properties of the system. 

However, in the case of molecules in solution, the environment is highly dynamic. 
The motion of the environment will lead to time dependence of the transition 
frequencies and couplings in the system (fluctuations), as well as to the exchange 
of energy with the environment (dissipation). The eigenstates of the system, found 
by diagonalizing the Hamiltonian for a given configuration of the bath, will only exist 
for a short time. As soon as the bath changes its state, the original eigenstates will 
mix. This situation can be called a dynamic bath. 

The effect of such a dynamic bath on the charge transport has been calculated 
from molecular dynamics simulations combined with quantum chemical methods 
[HI [T5] . Dynamic disorder (in the form of conformational fluctuations) is found to 
be important, as models with only static disorder cannot realistically describe charge 
transfer in DNA. These studies also strenghten the expectation that correlations in 
the fluctuations might be important in DNA. In the case of excitons, the coupling to a 
dynamic environment has been shown to be reflected in non-linear optical observables 

m- 

The interplay of electronic coupling with fluctuations and dissipation, induced 
by an environment, is regarded as an important problem in chemical physics. The 
resulting energy transport is traditionally modelled with Forster theory, but this 
cannot describe the energy transfer in the intermediate coupling limit or the presence 
of correlated fluctuations [17] . In the context of quantum networks [18] , it has been 
realized that the presence of fluctuations can increase the efficiency of excitation 
transport [191 [lOl ■ Recently, using two-dimensional optical spectroscopy, it has been 
discovered that quantum coherence in photosynthetic light-harvesting complexes can 
live much longer than expected from perturbative treatments of the system-bath 
interaction [21]. Although it was previously suggested that the origin of this effect 
might be found in correlated disorder, more recent studies have found that long-lived 
coherence can originate in a proper treatment of the time scale of bath fluctuations |22j . 
These treatments have, however, been limited to uncorrelated disorder. The effect of 
interplay between a dynamic bath and correlated fluctuations is yet unknown. 
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Figure 1. Structure of the six base pair DNA helix and of a single base pair, with 
the transition dipoles indicated. The DNA structure was generated with 3DNA 

This highlights the importance of a proper understanding of the quantum 
dynamics of a system in contact with its environment, including the possibility of 
correlated fluctuations. The effects of an environment that evolves on a time scale 
comparable to the excitonic dynamics, as well as the presence of environmental modes 
that couple to different chromophores simultaneously, remains to be studied. Rigorous 
theories have been developed to deal with the dynamics of a system in contact with a 
dynamic, quantum mechanical bath. The hierarchy of equations of motion approach 

[Ml US [13 m im HH] can be used to include dynamical fluctuations, as well 
as dissipation. Numerical path integral approaches that can deal with the dynamic 
effects induced by a quantum mechanical bath were introduced [29l [30l [311 [32] , and 
applied to the charge transfer in pieces of DNA [S^ . 

In this paper, we investigate the exciton dynamics in DNA in the presence of 
a dynamic bath. We focus on the effect of correlations in the fluctuations on the 
dynamics, and show how their effects are observable in the non-linear optical response. 
In particular, we calculate the anisotropy decay as a function of a coherence time and a 
population time, which has been used in experiment to characterize exciton transfer in 
conjugated polymers |34| . We introduce the theory, which includes a proper treatment 
of the environmental time scale as well as the presence of correlated fluctuations in 
sectional In section [3] we present the resulting exciton dynamics and calculated linear 
absorption spectra and two-time anisotropy decay data. Section [J] contains concluding 
remarks. 

2. Model 

2.1. Hamiltonian and equations of motion 

In this paper, we consider a piece of six base pairs of poly(dA) poly(dT) B-DNA, 
which has been studied experimentally and theoretically as a model compound. In 
this structure, shown in flgure[Tl each base pair contains a dA and a dT base, and all 
base pairs lie flat in planes perpendicular to the helix axis. The relative orientation of 
two adjacent base pairs is deflned by the twist angle. Although this angle can fluctuate 
considerably in solution, our model assumes a flxed twist angle of 36 degrees, which 
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is the average value for B-DNA. For this value, the helix makes a full turn every 10 
base pairs. 

In each base, we include a single strong ultraviolet transition, which is modelled 
as a two-level system. The frequencies and transition dipoles are chosen following |35| . 
The transition in the dT base is between the 5*0 and states, and has a frequency of 
er = 37.500 cm~^. In the dA base, the strongest absorption is found for the S'o — >■ S2 
transition, with a frequency of ea = 38.800 cm~^. Both transitions have a transition 
moment of 3. 71? in water. The transition dipoles lie in the plane of the bases, but 
their directions are not well known. Experimental and calculated results vary over 
about 90 degrees for adenine, and 30 degrees for thymine. To obtain a definite model, 
we adopt the experimentally determined values, as shown in figure [T] In adenine, the 
dipole flnA makes an angle of 66 degrees with the C-NH2 bond. In thymine, we use 
+14 degrees from the C=0 bond for the dipole jlnT- This results in an angle of 117 
degrees between the two transition dipoles in each base pair. 

The base pairs are labelled with indices n and m. The standard exciton 
Hamiltonian, which contains the excitation energies of the 2A'' bases, as well as 
interactions between them, is given in terms of the creation and annihilation operators 
and c by 

JV JV 

Hs = £iC„jC„i + Jni,mjCj^iCmj- (1) 

n=lie{A,T} n.m=li^je{A,T} 

In the B-DNA structure, the couplings were estimated to be JnA,nT = 248 cm~^, 
JnA,(n+i)A — 217 cm~^ and JnT,{n+i)T = 170 cm~^ [3S]. The couplings between 
stacked T bases and between the A and T in the same base pair agree well with 
results from quantum chemical calculations. Using these methods, however, a much 
larger value (872 cm"-'^) was found for the interaction between stacked A bases |36) . 
This 5 times larger value would significantly increase the length over which coherence 
can be present in the stack of A bases. In addition, the effect of interactions between 
non-adjacent bases must be considered. Here, we will not try to model the coupling 
strengths in more detail, but restrict our discussion to the parameters given above. 

In addition to the direct interactions between the base pairs, the excitations 
in DNA strongly interact with the environment, as can be seen from the Stokes 
shift observed in experiment [B]. The environment is modelled as a collection of 
independent harmonic oscillators, labelled by an index a, with masses m^, coordinates 
Xon momenta pa and frequencies Wq,, which gives the bath Hamiltonian Hb — 
Set (-Pq/2"^c( + TTia^a^a/'^) ■ The interaction between the bases and the environment 
is given by 

JV 

HSB=~'^ ^ '^gm.aVniXa, (2) 
n=li6{A,T} a 

where Vni = cl^^Cni- This system-bath interaction leads to fluctuations in the excitation 
frequencies of individual bases, and to dissipation of energy into and out of the bath. It 
describes the dynamics of excitons in the system, including their relaxation within the 
excitonic bands. We do not include couplings of the bath to a single system creation 
or annihilation operator, which changes the number of excitons in the system, and 
eventually returns the system to the ground state. The information on the system 
bath coupling is contained in the parameters gni,a, which denote the effect of the bath 
mode with frequency tUa on the transition energy of the ni"^ base. Because the bath 
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modes are harmonic oscillators, all information on the coupling of the system to the 
bath is contained in the spectral densities 

Jn^.raM = ^ ^^^^^^<5(c. - C.„). (3) 

— ' ZrUaUJa 

a 

A useful model for the spectral densities, which employs a single bath time scale 
for each base, is given by the overdamped Brownian oscillator. 

In principle, our treatment is not limited to overdamped modes, the approach used 
here can be extended to treat a more general Brownian spectral distribution [37l [38] . 
In the overdamped case, and assuming the high temperature limit, the correlation 
functions for the effective bath modes are given by [23^ .25] 

T .(f\ — „ . .p-lr^i.mM (K\ 

with 

2 

The functions Lni,ni{t) are autocorrelation functions. Their real parts describe the 
magnitude and the time scale of the fluctuations in the ni*^ base energy, determined 
at a given temperature T = l/ksP by \ni,ni and ^ni,ni, respectively. The imaginary 
parts of the correlation functions are responsible for the dissipation of energy. The 
presence of correlations between fluctuations on different bases is modelled by the 
cross correlation functions Lni^mj for 7^ f^j- They determine the degree to 
which fluctuations in the ni*'^ and the mj**^ frequency are correlated. If these 
quantities are all zero, all base energies fluctuate independently. In the other extreme 
case, the fluctuations are perfectly correlated if Xni,mj = Xni,ni = Xmj.mj and 
lni,mj = lni,ni = '~fmj,mj- In this casc, the ui^^ and mj*'^ base are coupled to a 
common bath, with the same strength. 

The dynamics generated by the complete Hamiltonian, which we will call the 
Hamiltonian of the medium, Hm = Hg + Hb + -ffsB, is given by the Liouville equation 

R{t) - ~iHl,R{t), (7) 

where B = [A,B] and R{t) is the density matrix that contains all the degrees of 
freedom of the system as well as the bath. Note that at this point the equation of 
motion is entirely equivalent to the Schrodinger equation, where the same Hamiltonian 
Hm describes the time evolution of a wave function. The formal solution of ([7]) is given 
in terms of the propagator G{t;to) = exp{—iH^,j{t — to)) by R{t) = G{t;to)R{to). 
This allows, in principal, the calculation of the density matrix R{t) at all times. The 
interesting observables are, however, normally defined only in terms of the 2N degrees 
of freedom of the system, and do not require the knowledge of the infinite-dimensional 
density matrix R{t). One therefore introduces the reduced density matrix, which 
only includes the system's degrees of freedom, as p{t) = tiBR{t) — J2ai^\^('^)\'^) ■ 
The time evolution of this reduced density matrix for a system linearly coupled to a 
harmonic bath can be found using numerical path integral techniques [SUl [SB [2S] • An 
alternative and efficient method, which we will use here, employs a set of equations of 
motion for the system's reduced density matrix and multiple auxiliary density matrices 
m m [13 [21 [H [m [2H1 [3H]- it can be derived starting from the Feynman- Vernon 
influence functional approach [23l [25] . The reduced density matrix for the electronic 
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transitions is given by path integrals over coordinates (j){t) and (j)'(t). The influence 
functional at time t, when the system was prepared in a factorized state at time 0, 
can be calculated by generalizing the procedure described in [55], and is found to be 

Fmt)}] = E f'^t'f dt"K><[{0(t')}](ReL™,™,(t'-t")V^„^, [{</>(*")}] 

ni,riij ^ ^ 



2lmL^,,^j{t' -t")V,:,mn}]). (8) 



where {(/>} denotes the pair (j) and 0', Fni['?^(0] the representation of the operator 
Vm, V;^[{<^(t)}] = Vmim - Vnd^'it)] and V°d{<f>it)}] = VmlHt)] + VmWit)]. 

Assuming commuting operators Vni, a hierarchy of equations of motion can be 
derived by introducing auxiliary density matrices, indexed by a set of indices Uss' ■ 
The reduced density matrix for the electronic transitions is found by setting all the 
indices to zero. It nonperturbatively contains the effects of the bath, as ensured by 
the presence of the auxiliary density matrices. In the high-temperature approximation 
[23], the hierarchy including cross-correlation terms, is given by 



i E Uss' (< 



ss' / 

c.,.K.p":.' it) - cL,p":=' {t)Vs') - z E v;v<- 



ss ss 

where we have introduced the notation s and s' for the pair ni and nf^, — Ugs' ±1- The 
coefficients Css> are given by ([6]). For systems at lower temperature, the hierarchy can 
be extended with low temperature correction terms [24i 125) . The lowest member of 
the hierarchy, which corresponds to the physical reduced density matrix, fully includes 
the coupling of the system to the bath. The deeper layers can be understood as 
bookkeeping devices which store the state of the bath, and, importantly, the coherences 
between system and bath states, at earlier times. 

To simplify the physical picture as well as the numerics, we will discuss only 
cases in which the fluctuations on two bases are either uncorrelated, or completely 
correlated. We will furthermore assume that the strength and the time scale of the 
fluctuations in each base are the same. In this case, the equations of motion simplify 
to the form 



pi-Ht) 



: E ns {cV.p"^ - c*p-^ Vs)-zY,K P"^ , (10) 



with c = A(— 17-1- ^). The sum over s runs over the terms in the system bath coupling. 

If all fluctuations are uncorrelated, s denotes a base ni, and Vs = c\^Cni- The model 
then includes 12 system bath operators V for a six base pair helix. In the case where 
the fluctuations in tA and ct on the same base pair are perfectly correlated, but 
uncorrelated with the frequencies on other base pairs, s refers to a base pair n, and 
Vs = cl^j^CnA + cliT'^riT- Wc furthermore approximately include the low temperature 
terms as described in [40 . The dynamics obtained from this equation of motion fully 
includes the fluctuations and dissipation induced by the interaction with the bath, 
without relying on a perturbative or fast bath approximation. 
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In the simulations presented in this paper, we use A = 2.5 • 10'^ cni^^ and a bath 
time scale tb — I/7 = 50 fs. The choice of A, which is an order of magnitude larger 
than the couplings J, is motivated by the observed large Stokes shift [SI The 
temperature is set to 300 K. For these values, jih'y = 0.5, which confirms the validity 
of the high-temperature approximation. 

2.2. Optical response 

2.2.1. Linear response. The dynamics of a single excitation in the DNA can be 
obtained directly by propagating ^ from a given initial condition. We now turn to the 
calculation of optical observables, which is formulated by combining the propagation 
the equation of motion with the correct sequence of matter field interactions. These 
are deduced by coupling the system to an external electric field E{t). The total 
Hamiltonian is the sum of the Hamiltonian Hm (defined as Hyi — Ha, + + i?sB)j 
the Hamiltonian for the field and an interaction term, which is given in the semi- 
classical description and in the dipole approximation by 

Hml = ■ E{t) = - ^ • E{t){ci + Cn). (11) 

n 

Note that in this equation we have slightly condensed our notation and used indices 
n and m to denote a single base rather than a base pair. This convention will be 
used throughout this section. The linear response requires the propagation of the free 
dynamics of the medium (i.e. in the absence of the field) after a single interaction 
with the light at a time to. The resulting polarization at time t can be expressed as 
the convolution of the electric field with a response function [41] . Assuming that the 
excitation is created by interacting with an ultra short laser pulse, the convolution 
becomes trivial, and the measurement of the polarization directly gives the response 
function as S{t) = {tr fip{t))o. The notation (• • ■)o indicates that the response 
function is averaged over all orientations of the sample with respect to the fixed 
coordinate system of the laboratory. The reduced density matrix is found from the 
perturbation expansion p(t) = i ti b G(t; to) ji^ R{tQ) (for t > io)- After making the 
rotating wave approximation (RWA), which is required for a consistent treatment of 
ultra short pulse excitation [41], it reduces to 

p{t)=iti-BG{t;to)fiR{to). (12) 

As can be seen, this reduces the number of terms in the linear response from two 
{jj,^ R = fiR — Rfi) to only one. 

Finally, the linear absorption spectrum is given as the Fourier transform of the 
response function (using t^ =0), 

1 1'°° 

A{u) = -Re / Ate'^^'fln ■ PmGno-Mt; 0). (13) 
-J Jo 

The factor 1/3 originates in the average over orientations. The indices n and m, which 
label individual bases, are understood to be summed over, and denotes the ground 
state. We have assumed that the system is in the ground state before interacting 
with the electric field. This is entirely reasonable in the case of optical excitation of 
DNA. The excitation energy is orders of magnitudes larger than the thermal energy, 
preventing thermal population of other states than the ground state. The interaction 
with the field then puts the system into a coherent superposition of a singly excited 
state and the ground state. Because the medium Hamiltonian conserves the number 
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of excitations, only the time evolution of such coherent superpositions is required for 
the calculation of the linear response. This is indicated by the notation G„o;mo(^;0)> 
where n and m denote the excitation of single bases, and 

(n|p(i)|0) =itrsG„O;™o(i;0)Mm-R(io). (14) 

It is essential to note that the propagator G does not only include the degrees of 
freedom of the system. It is written in the product basis of the system degrees of 
freedom and all bath states, and explicitly mixes the system and the bath through the 
action of Hsb- All mixing processes are included nonperturbatively by propagating 
the hierarchy of equations of motion starting from /l„i?(0), which boils down to the 
calculation of G„o;mo(^; 0)- 

2.2.2. Third-order response. Because the linear absorption contains no direct 
information about the exciton dynamics, and the second-order response vanishes in 
an isotropic medium, we will furthermore employ the third-order response |411 142) . 
Although it can be derived following the same procedure as outlined above for 
the linear response, the calculations is naturally somewhat more involved. More 
importantly, a third-order response function depends on three time intervals, and 
can be written generally as S{ti,t2,t3). A proper description of the system-bath 
interaction, which is not limited to an ultra fast bath, includes the memory stored in 
the bath of the system state at earlier times. The presence of this memory makes 
it impossible to split the reduced response functions, which depend only on the 
system's degrees of freedom, into functions that depend only on a single time variable, 
S{ti,t2,t3) ^ S{ti)S{t2)S{t3). The memory in the bath strongly affects the lineshape 
in nonlinear observables. Although it is quite naturally included in the simulation from 
a time-dependent Schrodinger equation, which can be used to calculate the nonlinear 
vibrational response (43l [44] , this approach does not describe the quantum mechanical 
nature of the bath. The approximation of a fast bath (often confusingly referred to as 
Markovian), commonly used in the case of electronic excitations neglects, as mentioned 
before, the bath memory, and therefore does not describe the partially inhomogeneous 
broadening observed in experiments, as has been shown in the case of a single spin 
system [45]. While the memory can be correctly included in the case of mostly static 
fluctuations [46], we shall develop here the full calculation for multiple coupled two- 
level systems in the presence of a bath which evolves on an arbitrary time scale. 

The third-order response can be measured in most detail by exciting the system 
with three short laser pulses, which interact at times ti, T2 and T3. In-between the 
pulses, the time evolution is dictated by the Hamiltonian Hm- We define variables for 
the evolution times by ti = T2 — ti, ^2 — T3 ~ ''"2 and t^ — t — t^. In the same way 
as in the linear response, the use of short excitation pulses allows the measurement of 
the third-order response function, which is given by 

^''^^1,^2,^3) = »'(trMG(i;r3)A^G(r3;T2)A^G(r2;Ti)A^i?(to))o-(15) 

Because the sample size is typically much larger than the wavelength of the light, 
the signal emitted in a third-order experiment is found in specific directions, given 
by linear combinations of the wavevectors of the incident pulses [17] . This allows the 
selection of a part of the response function. As mentioned before, in the impulsive 
excitation limit, it is necessary to apply the RWA for the system laser interaction. In 
the photon echo geometry, where the wave vector of the signal is equal to the sum of 
the wave vectors of pulse two and three, minus the wave vector of pulse one, three 
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contributions survive the RWA. The response function is the sum of these terms, which 
can be interpreted as a ground state bleach process, a stimulated emission process, 
and an induced absorption process. Their contributions to the response functions are 
given by 

S'seC^I, ^2, ^3) = i^il^m'" fin"Gm"'0;m"o{t: T'i) Pm" n") O i 

'5'lA(il, ^2, ^3) — ~ i^{Pw'n'''flwm''Gw'n'''-,wn''{t',T^)Pm''n'')o, (16) 

respectively, where we have defined 

(T = Pn' PnGon' ;0n{T2; Ti) (17) 

and 

Pm"n" — P7n'P7iGm''n''-,m'7i'{T3]T2)Gon'-,On{T2]Ti). (18) 

In these equations, repeated indices are understood to be summed over. The 
orientational average over the rank four tensor formed by the product of four dipoles 
can be written in Cartesian components as {jljljljT)o = ^ap-yS ^aiS-ySP^ tJ-^ tJ'^ IJ'"' ■ The 
tensor A contains the effects of polarization of the laser pulses. Three fundamentally 
different combinations of polarizations can be used, giving different values for the 
elements in A [35]. Combined with the sum over Cartesian components, labelled with 
a, fS, 7, S, they describe the averaging of the signal over an isotropic sample. 

During the coherence time ^3, the induced absorption contribution contains 
coherent superpositions between one- and two-exciton states. These two-exciton states 
are labelled with indices w and w' in the last line in (IT51) . We work in the site basis, 
where the two-exciton basis states are given as the direct product of two excitations. 
For example, a two-exciton state can be written as \w) = \w1W2), in which case base 
Wi and W2 are excited, while the electrons in all other bases are in the ground state. 
Consequently, the transition dipoles between one- and two-exciton states are given by 

t^wn — P'wlw2^n ^wl.7iPw2 ^" ^w2,npwl- 

It is important to note that the above-mentioned memory stored in the bath 
is fully included in these expressions for the nonlinear response. The propagators 
G are matrices in the product basis of system and bath degrees of freedom, and 
multiplications over the bath degrees of freedom are understood by writing the explicit 
dependence of the propagators on two times. In a more complete but complicated 
notation, one might write for example 

Pm"f" = A^m'A^nG,^'^i'';m'„'(^2)Gj„^0ri(^l)-^0 ' (l^) 

where S, T and U denote basis states in the Liouville space of the bath degrees of 
freedom. S and T are understood to be summed. This dependence of the propagators, 
which describes the memory of the bath and of the system-bath coherence, is correctly 
included by propagating the hierarchy of equations of motion. The simulation 
procedure of, for example, the stimulated emission contribution, is as follows. First, 
the matrix elements of the reduced density matrix, as well as all the auxiliary density 
matrices, between a singly excited state and the ground state, are populated according 
to their transition dipole. The hierarchy of equations of motion is then solved 
for a time ti. The second multiplication with the transition dipoles gives matrix 
elements of all the (auxiliary as well as physical) density matrices between two singly 
excited states. The hierarchy for these matrices is propagated for a time <2- A third 
multiplication with the dipoles again gives a coherence between singly excited states 



Correlated fluctuations in the exciton dynamics and spectroscopy of DNA 



10 



and the ground state, which is propagated for a time t^, followed by a final transition 
dipole multiplication and a trace operation to obtain the signal. This procedure is 
repeated for all required combinations of Cartesian components. 

While the linear absorption probes only the dynamics of coherent superpositions 
between the ground state and one-exciton states, the third-order response function 
contains the dynamics of matrix elements {n\p\m) during t2- These terms describe 
the motion of a single excitation in the system, as well as the time evolution of quantum 
coherences between system states. The excitation dynamics can therefore be measured 
as a function of ^2, and be correlated with the time evolution of the coherences during 
t\ and tz- 



2.2.3. Two-time anisotropy A useful measurement of this information in a system 
where the transition dipoles of individual bases point in various directions, as is the 
case in DNA, is the two-time anisotropy decay. It can be found by setting ^3 = 0. 
In this case, the stimulated emission is cancelled by a contribution from the induced 
absorption, and the signal simplifies considerably. 



S{tl,t2) = ( p.mP-m<y — ^ Ppi^pPn 



(20) 

pjtn I Q 

This anisotropy can be measured independently in the parallel polarization geometry 
(using four pulses all polarized in the same direction), which gives a signal S'P'*'', and 
in the perpendicular geometry (where the final pulse pair is polarized perpendicularly 
to the first pair), which yields In simuations, these two reponse functions are 

found by choosing the appropriate coefficients A. The two-time anisotropy is then 
calculated as 

Q U , ^_ gP"(^l,^2)-^P"(tl,t2) 

The anisotropy describes the average rotation of the dipole of the excitation during 
the times ti and t2- It is a complex valued quantity, and in the rest of this paper we 
will only consider its absolute value. In numerical calculations, it is an advantage that 
the propagation of the two-exciton states is not required. 



2.3. Static limit 

The time scale of the bath tb is modelled by the parameter 7 = I/tb- The static 
limit is recovered by assuming an ultra slow bath, characterized by the limit 7 — >■ 0. 
The correlation function becomes 

m = J. (22) 

We see that the imaginary part of the correlation function vanishes in this limit, while 
the real part is a constant. In physical terms, this limit means that fluctuations in 
the base transition energies are now frozen. The system can be viewed as an ensemble 
of helices, in each of which the bases have different energies. In addition, in the case 
of incomplete correlation, each base within a given helix has a different transition 
frequency. The system is therefore no longer ergodic, as should be expected in the 
static limit. The absence of an imaginary term in the corrcilation function shows that 
there is no dissipation. This can be explained by considering that the bath modes are 
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now ultra slow; they have too low frequencies to accept any energy from the system 
in a finite time. 

In the static limit, the linear absorption spectrum can be calculated from Fermi's 
golden rule, 



where E]^ is the energy of the k eigenstate with eigenfunction (f>k, and fk = 
(Sn '/'fcn/Jn)^ IS its osciUator strength. The average, indicated by the notation (• • •) 
is over realizations of the static disorder. Similar expressions can be derived for the 
non-linear response. 

The dynamics of a single excitation in the site basis can be found by projecting 
the initial condition on the eigenstates, propagating the density matrix in the 
eigenbasis and projecting back. In each realization, an initial population on site n 
gives a density matrix Pkq{0) — 4>kn4>qn- Its time evolution is given by Pkq{t) = 
exp(—i(Ek — Eg)t)pkg{0), and the final density matrix in the site basis is obtained as 
Pnmit) = J2kq'^kn4'qmPkqi't)- This density matrix can be averaged over realizations 
to obtain the average population on a given site, or the coherence between sites. 

3. Results and Discussion 

3.1. Delocalization 

The properties of elementary excitations in the DNA are determined by a competition 
between electronic interactions between bases, which favour extended states, and 
environment-induced fluctuations, which lead to localization. In the case of DNA, 
the coupling to the environment, which is determined by the reorganization energy 
and the temperature according to (|22|) . is an order of magnitude larger than the 
electronic couplings. We therefore analyse our results in the site basis, and discuss 
populations on individual bases and coherences between bases. One might wonder if 
the large reorganization energy does not prohibit any significant quantum coherence 
between bases. In the static limit, the properties of the system can be calculated 
by averaging over all members in the ensemble, or, equivalently, over individual 
realizations of the disorder. For each realization, the Hamiltonian can be diagonalized 
to yield the eigenstates. These eigenstates are used to calculate the observable of 
interest, which is finally averaged over the realizations of the disorder. The localization 
length in a linear chain can be estimated using the properties of the eigenstates 
[49j [50] . Such an estimate would lead to the conclusion that excitations are spread 
over a few bases. Here, we directly calculate the coherence between two bases in an 
ensemble of dimers, each consisting of only two neighbouring A (or, equivalently, T) 
bases. For a particular system in the ensemble, the Hamiltonian can be written as 
H = (<5e/2)(4c2 - c[ci) + J{c[c2 + cjci). 

While the value of the coupling J is the same throughout the ensemble, the energy 
difference between the two bases 6e varies from one system to another. This difference, 
which originates in the fluctuations, is characterized with a Gaussian distribution, 
which has a standard deviation cr ~ y/2X/ (3. In each realization, diagonalization of 
the Hamiltonian gives the eigenstates \+) and |— ) (note that these are the symmetric 
and anti-symmetric states only in the absence of disorder). For a positive value of (5e, 
the state with the highest energy is found mostly on base 2, while for negative de it 




(23) 



k 
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Figure 2. Linear absorption spectra A(ijj) in (a) the static limit and (b) with 
dynamic fluctuations for uncorrelated baths (solid line), correlated fluctuations in 
each base pair (dotted line) and fully correlated fluctuations (dashed line). 



occupies mainly base 1. In either case, the coherence C — (l|+)(+|2) is weU-defined, 
and given by 



C^ihhl^^ (24) 
D-SeVD 



with the energy gap between the eigenstates D = %/ Se'^ + 4J2. Integrating over a 
Gaussian distribution, we find the average coherence in the presence of disorder to be 

iO = :lIfeW^y'Ko{{J/an (25) 

V ZTT 

where Kq is a modified Bessel function of the second kind. In the absence of disorder, 
which amounts to taking the hmit of a/ J to zero, the coherence reduces to the value 
for completely delocalized states, (C) = 0.5. In the other limit, where a/ J — > cxi, 
complete localization destroys the coherence, which approaches zero. Surprisingly, 
however, substantial coherence is present for large values of the disorder. In particular, 
for a/ J = 10, we find (C) = 0.19. This shows that, even though the coupling to the 
environment is quite strong compared to the electronic interactions between the bases, 
the excitations cannot be understood as completely localized. Even though, based on 
the estimate mentioned before as well as on earlier work the localization length is 
not expected to be larger than a few bases, the coherence between bases cannot be 
ignored in a description of energy transport. 



3.2. Linear absorption spectrum 

The linear absorption spectrum, calculated from (I13p . is shown in figure[2]for the static 
as well as the dynamic environment. The static results were calculated by averaging 
over 10^ disorder realizations. We observe that the line shape appears similar in all 
three models, but that the correlations lead to a shift in the absorption maximum in 
the case of a static environment. This can be understood from the decrease of the 
localization in the presence of correlations. For completely correlated disorder, the 
eigenstates are delocalized over all available bases. The delocalization changes the 
character of the bright states, which leads to a shift in the absorption spectrum. 



Correlated fluctuations in the exciton dynamics and spectroscopy of DNA 13 
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time (fs) time (fs) 

Figure 3. Population on tiie first A base, p-iAiAi^)^ after initial excitation of 
this base for uncorrelated fluctuations (solid line), correlated fluctuations in each 
base pair (dotted line) and completely correlated fluctuations (dashed line). 

The spectra for a dynamic bath are shown m panel (b). The spectra for correlated 
and uncorrelated fluctuations appear almost completely similar; the shifts observed in 
the static spectrum disappear. We conclude that the presence of correlation cannot 
conclusively be established from the linear absorption spectrum. 

3.3. Exciton dynamics 

Although the difference between correlated and uncorrelated fluctuations can't be seen 
in the linear absorption, it has a strong effect on the exciton dynamics. This can be 
clarified by performing a numerical experiment, in which the first A base is artificially 
excited. The subsequent time dependence of the population on this site is shown 
in figure [31 The fully correlated case is rather trivial. In this case, the fiuctuations 
have no effect at all, and the dynamics are generated by the system's Hamiltonian 
without influence from the bath. The time scale of the bath is consequently not 
important, the time evolution is the same for all values of 7. For perfectly correlated 
fiuctuations in the A and T base in each base pair, but no correlations between the 
base pairs, the exciton dynamics is, perhaps surprisingly, slower than for completely 
uncorrelated fluctuations. This effect is only present in the dynamical model for the 
environment. As expected, the dynamic and static baths give rise to almost the same 
exciton dynamics in the flrst 100 fs, while significant differences appear for longer 
times. This is especially clear in the population dynamics in the case of uncorrelated 
fluctuations. 

The coherence between two neighbouring A bases is plotted in flgure [H The 
coherence is clearly increased by the presence of correlations along the stack. The 
real part of the coherence (not shown) is much larger in the presence of dynamic 
fluctuations as compared to a static bath, when it is zero for a dimer. We flnd similar 
results for the coherence between the A and T bases (not shown). 

3.4^. Two-time anisotropy decay 

The time dependence of the anisotropy, as measured in a third-order non-linear 
experiment, contains information about the exciton dynamics, and has been used 
for that reason to study electronic excitations in conjugated polymers [34], as well as 
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Figure 4. Absolute value of the coherence between neighbouring A bases, 
PiA2A(t)i after initial excitation of the first A for uncorrelated fluctuations 
(solid line), correlated fluctuations in each base pair (dotted line) and completely 
correlated fluctuations (dashed line). 



vibrational modes in peptides [ST] . The presence of disorder has been shown to lead to 
relaxation of the anisotropy within 100 fs in conjugated polymers [S^. This effect could 
be explained by interaction with an environment which contains either a slow (static) 
or a fast mode, compared to the time scale of the system's dynamics as probed in the 
experiment. In the presence of static disorder, interaction with light creates excitation 
of eigenstates of the system, as well as coherent superpositions of these eigenstates. 
The populations in the eigenbasis do not evolve under the system Hamiltonian, while 
the coherences pkq only collect a phase factor determined by the energy difference 
Ek — Eg. Incoherent relaxation is introduced by the fast environmental modes, which 
are treated perturbatively. Because of the assumption that these modes relax much 
faster than any time scale in the system, a master equation can be derived for the 
populations of each eigenstate. Such a description has been successfully applied to 
understand the dynamics of excitons in self-assembled aggregates at a temperature 
low enough to freeze the solvent, in which case the assumption of a static environment 
is valid ^551. 

If this type of calculation is applied to our current results of the anisotropy 
relaxation in the first 100 fs, this assumption would imply that the fast environmental 
modes have characteristic dynamics on time scales considerably smaller than 10 fs. 
Here, we rather model the decay of the anisotropy using a single time scale tb of the 
bath. For times smaller than tb, a static description of the dynamics is valid. The 
dynamics are then governed only by coherent evolution under the system Hamiltonian 
in each static realization. For times comparable to or larger than tb, the bath time 
scale leads to the destruction of the initial eigenstates, and the results will deviate 
from a picture with static disorder only. 

Because our model does not include the (presumably slow) rotation of the DNA 
molecule, the anisotropy does not decay for a single base. All the decay in the 
anisotropy can therefore be ascribed to energy transport in the DNA. The anisotropy 
decay in the DNA model is quite complex. It involves transfer between the A and T 
bases in a pair, which leads to a large (117 degrees) rotation of the transition dipole, 
as well as exciton migration along the stack, which is energetically more favourable, 
but only rotates the transition dipole by 36 degrees in a single hop. It is not a priori 
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Figure 5. Two-time anisotropy decay Sttad (il > ^2) for (a) a single base pair with 
uncorrelated fluctuations (b) a single base pair with correlated fluctuations, (c) 
two neighbouring A bases with uncorrelated fluctuations and (d) two neigbouring 
A bases with correlated fluctuations. 



clear which mechanism will contribute more to the anisotropy decay. To understand 
these two contributions separately, we present the decay of the anisotropy in a single 
base pair, as well as in two neighbouring A bases in figureO In general, the anisotropy 
decays faster for transfer between the A and T base, for both correlated as well as 
uncorrelated fluctuations. This indicates that in DNA, the effect of the larger dipole 
rotation is more important than the difference in the A and T excitation energies. 
In the case of fully correlated fluctuations, however, the recurrence time in the AT 
base pair is shorter than in the AA stacked pair, as can be seen in the cross sections 
in figure ID Furthermore, the presence of fluctuations always slows down the initial 
decay of the anisotropy, as can be seen from the faster decay in the case of correlated 
as compared to uncorrelated fluctuations. 

We now turn to the discussion of the two-time anisotropy decay in DNA, 
calculated with a bath time scale of 50 fs. Comparing the results with the decay of the 
anisotropy in a single base pair with uncorrelated fluctuations (for clarity repeated in 
figure [7] (a)), we see that the decay is faster for the extended helix. The excitation now 
has the freedom to move along the stack as well as between base pairs, which explains 
the faster decay. We next analyse the effect of correlations in the fluctuations, which 
clearly leads to differences in the two-time anisotropy results. 

For ti — 0, plotted in figure [S] (b), the fully correlated and the AT correlated 
cases show similar results. In this case, the transfer within a base pair is completely 
coherent, as can be seen from the presence of a minimum near t2 = 80 fs, after which 
the anisotropy increases again. This is the same coherent recurrence as observed in 
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Figure 6. Slices through the two-time anisotropy for a single base pair with 
uncorrclated fluctuations (solid line) and correlated fluctuations (dashed line) as 
well as for two neighbouring A bases with uncorrelated fluctuations (dash-dotted 
line) and correlated fluctuations (dotted line). 
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Figure 7. Two-time anisotropy decay SttadC^i i ^2) for (a) a single base pair, 
(b-d) 12 bases with (b) uncorrelated fluctuations in each base, (c) fully correlated 
fluctuations in each base pair, but no correlations between base pairs and (d) fully 
correlated fluctuations. 



the populations and coherences before. In the presence of uncorrelated fluctuations, 
such recurrences are not present. In this case, the anisotropy decays steadily, in line 
with the behaviour of the population in a dynamic bath shown in figure |3] (a) . We 
conclude that the anisotropy measurement contains clear information on the presence 
of correlated fluctuations. 



Correlated fluctuations in the exciton dynamics and spectroscopy of DNA 

0.40 



-| 1 r 




J I L 



150 
t, (fs) 



200 



17 



Figure 8. Slices through the two-time anisotropy for a single base pair (dash- 
dotted line) and a helix with 12 bases and uncorrelated fluctuations in each base 
(solid line), correlated fluctuations within each base pair (dotted line) and fully 
correlated fluctuations (dashed line). 



Because our simulations treat only a small piece of DNA, one might expect finite 
size effects to play a significant role. To address this, we have repeated our simulations 
for a helix with five base pairs, and found almost identical results. We therefore 
conclude that finite size effects do not significantly change our conclusion and believe 
that our study of the short time dynamics provides insight in the properties of real 
DNA samples. 



4. Concluding Remarks 

In conclusion, we have modelled the exciton dynamics and non-linear optical response 
in a model BDNA helix, in the presence of a dynamic bath. In this study, we have 
used the hierarchy of equations of motion approach, which nonperturbatively includes 
the environment-induced fluctuation and dissipation, and extended it for the first time 
to include correlations in the fluctuations. The correlations have only a small effect 
on the static linear absorption spectrum, and are not visible in the linear absorption 
in the presence of a dynamic environment. We find that the presence of correlation in 
the fluctuations does have a strong effect on the exciton dynamics, which is reflected 
in the two-time anisotropy decay. 

Our current simulations are limited to a small piece of the DNA helix. We have 
found that the system size does not significantly influence our results. However, for 
the treatment of other systems, an extension of the method presented here to a larger 
system size is desirable. Although this is in principle straightforward, the simulation of 
larger systems is computationally expensive. In the case of DNA, the high temperature 
approximation is applicable, which speeds up the computation. However, because of 
the large ratio of the reorganization energy to the coherent coupling, the number 
of tiers in the hierarchy must be rather large, especially if we want to extend the 
calculation to a slower bath. The deeper tiers include terms which involve the action 
of multiple system-bath operators on bases that are separated in space. Correlations 
between the excitations on bases which are separated by more than the localization 
length are expected to be small. Computational improvements based on this physical 
argument, be it automatic truncation of the hierarchy as proposed recently [54) . or 



Correlated fluctuations in the exciton dynamics and spectroscopy of DNA 



18 



new algorithms that take into account the locahzation of exciton states exphcitly, are 
desirable. Furthermore, charge transfer states may be important in DNA [71[3S], and 
could be included in further development of the calculations presented here. 

The helical structure is important for a correct description of the anisotropy decay. 
For a precise comparison with experimental results, it is therefore desirable to include 
variations in the structure. Finally, while this work was limited to poly(dA)poly(dT) 
DNA, a straightforward extension of the method to include other base sequences would 
help in the interpretation of experiments. 

References 

[1] F. Fiebig. J. Phys. Chem. B, 113:9348, 2009. 

[2] D. Onidas, T. Gustavsson, E. Lazzarotto, and D. Markovitsi. Phys. Chem. Chem. Phys., 9:5143, 
2007. 

[3] I. Buckvarov, Q. Wang, M. Raytchev, A. Trifonov, and T. Fiebig. Proc. Natl. Sci. USA, 
104:4794, 2007. 

[4] B. Bouvier, J. P. Dognon, R. Lavery, D. Markovitsi, P. Millie, D. Onidas, and K. Zakrzewska. 

J. Phys. Chem. B, 107:13512, 2003. 
[5] S. Tonzani and G. C. Schatz. J. Am. Chem. Soc., 130:7607, 2008. 

[6] D. Markovitsi, D. Onidas, F. Talbot, S. Marguot, T. Gustavsson, and E. Lazzarotto. J. 

Photochem. Photobiol. A, 183:1, 2006. 
[7] E. R. Bittner. J. Chem. Phys., 125:094909, 2006. 
[8] E. R. Bittner. J. Photochem. Photobiol. A, 190:328, 2007. 
[9] D. Porath, A. Bezryadin, S. de Vries, and C. Dekker. Nature, 403:635, 2000. 
[10] A. V. Malysliev, E. Diaz, F. Dommquez-Adame, and V. A. Malyshev. J. Phys.: Condens. 

Matter, 21:335105, 2009. 
[11] S. Roche. Phys. Rev. Lett., 91(10):108101, Sep 2003. 

[12] P. Carpena, P. Bernaola-Galvan, P. Ch. Ivanov, and H. E. Stanley. Nature, 418:955, 2002. 
[13] E. Diaz and F. Dommguez-Adame. Chem. Phys., 365:24, 2009. 

[14] R. Gutierrez, R. A. Caetano, B. P. Woiczikowski, T. Kubar, M. Elstner, and G. Cunibcrti. Phys. 

Rev. Lett., 102(20):208102, May 2009. 
[15] B. P. Woiczikowski, T. Kubaf, R. Gutierrez, R. A. Caetano, G. Cunibcrti, and M. Elstner. J. 

Chem. Phys., 130:215104, 2009. 
[16] K. Hyeon-Deuk, Y. Tanimura, and M. Clio. J. Chem. Phys., 128:135102, 2008. 
[17] D. Beljonne, C. Curutchet, G. D. Scholes, and R. J. Silbey J. Phys. Chem. B, 113:6583, 2009. 
[18] J. Cao and R. J. Silbey J. Phys. Chem. A, 113:13825, 2009. 

[19] M. Mohseni, P. Rebentrost, S. Lloyd, and A. Aspuru-Guzik. J. Chem. Phys., 129:174106, 2008. 
[20] F. Caruso, A. W. Chin, A. Datta, S. F. Huelga, and M. B. Plenio. J. Chem. Phys., 131:105106, 
2009. 

[21] G. S. Engel, T. R. Calhoun, E. L. Read, T. K. Ahn, T. Mancal, Y. C. Cheng, R. E. Blankenship, 

and G. R. Fleming. Nature, 446:782, 2007. 
[22] A. Ishizaki and G. Fleming. J. Chem. Phys., 130:234111, 2009. 
[23] Y. Tanimura and R. Kubo. J. Phys. Soc. Jpn., 58:101, 1989. 
[24] A. Ishizaki and Y. Tanimura. J. Phys. Soc. Jpn., 74:3131, 2005. 
[25] Y. Tanimura. J. Phys. Soc. Jpn., 75:082001, 2006. 
[26] A. Ishizaki and Y. Tanimura. J. Phys. Chem. A, 111:9269, 2007. 
[27] J. Strumpfer and K. Schulten. J. Chem. Phys., 131:225101, 2009. 
[28] L. Chen, R. Zheng, Q. Shi, and Y. Yan. J. Chem. Phys., 131:094502, 2009. 
[29] C. H. Mak and D. Chandler. Phys. Rev. A, 44:2352, 1991. 
[30] M. Topaler and N. Makri. J. Phys. Chem., 100:4430, 1996. 
[31] N. Makri. J. Phys. Chem. A, 102:4414, 1998. 
[32] E. Sim. J. Phys. Chem. B, 108:19093, 2004. 
[33] H. Kim and E. Sim. J. Phys. Chem. B, 112:2557, 2008. 
[34] E. CoUini and G. D. Scholes. Science, 323:369, 2009. 

[35] B. Bouvier, T. Gustavsson, D. Markovitsi, and P. Millie. Chem. Phys., 275:75, 2002. 
[36] A. Czader and E. R. Bittner. J. Chem. Phys., 128:035101, 2008. 
[37] Y. Tanimura and S. Mukamel. J. Phys. Soc. Jpn., 63:66, 1994. 
[38] M. Tanaka and Y. Tanimura. J. Phys. Soc. Jpn., 78:073802, 2009. 



Correlated fluctuations in the exciton dynamics and spectroscopy of DNA 



19 



[39] A. J. Leggett, S. Chakravarty, A. T. Dorsey, M. P. A. Fisher, A. Garg, and W. Zwerger. Rev. 

Mod. Phys., 59:1, 1987. 

[40] A. Ishizaki and G. R. Fleming. Proc. Natl. Acad. Sci. USA, 106:17255, 2009. 

[41] S. Mukamel. Principles of Nonlinear Optical Spectroscopy. Oxford University Press, Oxford, 

1995. 

[42] M. Cho. Chem. Rev., 108:1331, 2008. 

[43] H. Torii. .J. Phys. Chem. A, 110:4822, 2006. 

[41] T.l.C. Janscn and J. Knocstcr. Acc. Chem. Res., 42:1405, 2009. 

[45] A. Ishizaki and Y. Tanimura. Chem. Phys., 347:185, 2008. 

[46] A. G. Dijkstra, T.l.C. Janscn, and J. Knocstcr. J. Chem. Phys., 128:164511, 2008. 

[47] S. Mukamel. Annu. Rev. Phys. Chem., 51:691, 2000. 

[48] R. M. Hochstrasser. Chem. Phys., 266:273, 2001. 

[49] V. A. Malyshev. J. Luminescence, 55:225, 1993. 

[50] A. V. Malyshev and V. A. Malyshev. Phys. Rev. B, 63:195111, 2001. 

[51] T. 1. C. Jansen and J. Knoester. J. Phys. Chem. B, 110:22910, 2006. 

[52] T. E. Dykstra, E. Hennebicq, D. Beljonne, J. Gierschner, G. Claudio, E. R. Bittner, J. Knoester, 

and G. D. Scholes. J. Phys. Chem. B, 113:656, 2009. 

[53] D. J. Heijs, V. A. Malyshev, and J. Knoester. Phys. Rev. Lett., 95:177402, 2005. 

[54] Q. Shi, L. Chen, G. Nan, R. X. Xu, and Y. Yan. J. Chem. Phys., 130:084105, 2009. 

[55] E. M. Conwell, P. M. McLauglin, and S. M. Bloch. J. Phys. Chem. B, 112:2268, 2008. 

[56] X. Lu and W. K. Olson. Nucleic Acids Res., 31:5108, 2003. 



